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During  a magnetospheric  substorm  the  plasma  in  the  E-region  of  the 
auroral  ionosphere  is  observed  to  be  unstable  (Booker,  1960).  The  instability 
is  due  to  a relative  drift  between  ions  and  electrons  in  a direction  per- 
pendicular to  the  geomagnetic  field  lines.  Such  a drift  can  occur  frequently 
in  the  auroral  ionosphere  at  an  altitude  of  approximately  110  km  because 
of  the  presence  of  an  electric  field  directed  normal  to  the  ambient  magnetic 
field.  In  this  region  the  motion  of  the  ions  is  somewhat  restricted  by 
collisions  with  the  neutral  particles,  while  at  the  same  time  the  lighter 
electrons  are  more  or  less  free  to  execute  a drifting  motion  with  a velocity 
equal  to  E x B/B^. 

Recently  Lee  e£  al.  (1971)  showed  that  a high  density  plasma  (N^  ^ 10^ 
cm“^)  is  susceptible  to  a high  frequency,  short  wavelength  type  instability, 
if  there  exists  a sufficiently  large  drift  velocity  between  ions  and  electrons 
perpendicular  to  the  background  magnetic  field,  called  the  high-frequency 
Hall  current  instability  (Lee  et  al . , 1971).  Their  linear  analysis  of  this 
Instability  results  in  a threshold  drift  velocity  V_  > 3(k_  T/M) , where  k„ 
is  Boltzmann's  constant,  T » T^  * T^  is  the  plasma  temperature,  and  M is  the 
ion  mass.  However,  observations  by  Kelley  and  Mozer  (1973)  seem  to  indicate 
that  nonlinear  effects  somehow  stabilize  the  waves  at  a speed  approximately 
equal  to  the  ion  acoustic  velocity  and  less  than  the  electron  drift  velocity. 

Following  the  work  of  Rogister  (1971)  on  the  equatorial  electrojet  we  give 
a nonlinear  analysis  of  the  high  frequency  Hall  current  two-stream  instability 
applicable  to  the  auroral  electrojet.  Our  results  show  that  the  unstable 


waves  indeed  travel  at  a speed  less  than  the  electron  drift  velocity.  Moreover 
we  made  an  order  of  magnitude  estimate  of  the  density  fluctuation,  and  this 
agrees  well  with  that  observed  by  Booker  (1960) . Also  the  proper  direction 
of  the  propagation  of  the  unstable  waves  is  shown. 


r 

2.  Discussion  of  che  Macroscopic  Equations 

In  che  region  of  che  auroral  electro jet  electron- ion  and  Ion-ion 
collisions  are  negligible  compared  with  collisions  between  neutrals  and 
charged  particles.  Moreover  ionization  and  recombination  processes  may 
be  ignored.  We  shall  use  che  two-fluid  equations  to  describe  che  simple 
model  of  che  electro jet  layer.  In  this  framework  the  fundamental  equations 
are 


If  + 7.(„  Vj,  . 0 


3V 


“"j  - j ’ ^ 


- 7p^  + nqj  (E  + 7 Vj  x B)- 


(1) 


(2) 


where  j ■ e,  i refers  to  the  electrons  and  ions  respectively,  and  Vj  is  Che 
collision  frequency  between  the  charged  particle  and  a neutral  particle. 
In  equations  (1)  and  (2)  charge  separation  and  finite  Larmor  radius  effects 
are  not  Included;  these  may  be  shown  to  be  negligible  in  che  region  of  Che 
auroral  electrojet. 

To  start  che  analysis  we  decompose  each  variable  according  to 


,j,  - 4,  + 5(>  (3) 

where  ^ is  che  spatial  average  of  and  iSi^  is  a fluctuation  component,  such 
that  its  average  <6<^>  *0.  In  accordance  with  the  Farley-Buneman  type 
instability  we  consider  electrostatic  waves  with  propagation  vector  per- 

/ 

pendicular  to  che  ambient  magnetic  field.  Also  we  shall  only  consider 
fluctuations  with  very  small  parallel  wave  numbers  as  these  are  more 
susceptible  to  resistive  instabilities  (Rosenbluth,  1965).  Furthermore, 
for  convenience  we  shall  assume  one-dimensional  propagation,  i.e.  we  let 
7 = §(3/3s),  where  ^ is  a unit  vector  along  the  direction  of  wave  propagation. 
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We  now  assume  Chat  the  average  quantities  are  in  equilibrium,  and  expand 
equation  (1)  as  follows: 


In.  0 

3c 


^ li'  ^ " -h 


Similarly,  expanding  eq.  (2)  one  obtains 


3V 


3t  * ^ Ti  «!!j'  * '’j% 


(n  + (Sn)m 


j 


+ -i-  (E  + -(V^  + 5V,)  + B] 
m^  - c -j  -J 


where  use  has  been  made  of  the  relations 


Taking  the  average  of  eq.  (2a)  and  noting  that 


<6Vj>  - 0 and  <(Vj  • §)-^  6V^>  - (V^ 


<6Vj>  - 0 


one  obtains 


3V 


S ^(5pJ 


+ <(3  • SV,)  5V  > •<-  ,,  > 

3t  -j  3s  -j  (n  + 6n)mj 


♦ t!  4 X 51  - V,  V, 


j -j 


Subtracting  eq.  (2b)  from  eq.  (2a)  one  obtains 


(4) 

(5) 


(2a) 


(2b) 


3 


36V. 


% > 


^ n,  - ^ «yj> 


(•  , X s' — +'^,-  j.  X \ ^ ~ (<5V  X B)  - V V. 

(n  + 6n)  nij  (n  + 6n)mj  c -j  - j “j 


Using  Cartesian  coordinates,  we  let  B ■ B(-l  ) and  § - s i + s i 

— 2 xxyy 

such  that  i ■ 0.  Talcing  the  vector  product  of  § with  eq.  (6)  and 

followed  by  the  scalar  product  with  i^,  we  obtain 


36V. 


n,  • § X + (V,  + 6V,)*  s 6V,  - <6V,  • S ~ 6V^> 


(6) 


-j  -i 


3s  -j 


-j 


3s  -j 


+ Vj  6Vj}  - s • 6Vj 


(7) 


q.B 


where  fl . = 


,th 


1.  = is  the  Larmor  frequency  of  the  j particle.  Also,  taking 

the  scalar  product  of  i with  eq.  (6),  and  summing  over  all  electrons  and 


ions,  we  obtain  the  following  equation 


)L.i  "j  ’ • ' 

■,L.i 


(8) 


where  we  have  used  the  simplification 


• <6V,  • S ^ 6V.>  - i <A  (g  . 6v,)2>  . 0 
-j  3s  -j  2 3s  -j 


From  the  Birkeland  current  system  we  note  that  the  field-aligned 
current  feeds  into  the  south  edge  of  the  auroral  arc  (see  Fig.  1). 
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y(NORTH) 


x(EASn 


From  the  current  continuity  equation  we  then  have 


J - 0 


iy  the  current  density  in  the  South-North  direction,  and  j,|  is  the 
Birlceland  current  which  feeds  into  the  auroral  arc.  The  current  in  the  y 
direction  is  largely  carried  by  the  electrons  on  account  of  their  high 
mobility. 


Now, 


(11) 


ne  V,  - ne  V 
iy  ey 


and  using  (3)  to  expand  (11),  and  taking  the  average,  we  obtain 


j - e(nV,  + <6n6V,  > - iiv  - <6n6V  >)  (12) 

•jy  ' ly  iy  Cy  £y 

since  the  current  In  this  region  Is  carried  mainly  by  the  electrons 
(V^y  <<  V^y) . It  follows  that 

1 a - neO 
Jy  ey 

Hence,  from  eq.  (12)  we  then  obtain 

<6n6V  > * n V.  + <6n6V.  > (13) 

ey  ly  ly 

We  are  now  considering  the  two-dimensional  flow  across  the  ambient  magnetic 
field.  With  the  magnetic  field  Bi^  and  the  primary  electric  field  in 
the  East-West  direction  given  we  have  a system  of  six  equations  in  the  seven 
average  fluid  variables  , n,  E^,  and  P.  Hence,  we  also  need  the  equation 
of  state  p » to  complete  the  description  of  the  system.  Here  it  has 

been  assumed  that  at  the  altitude  region  around  110  km  electrons  and  ions  are 
in  thermal  equilibrium,  i.e.  ■ T,  and  T is  considered  to  be  known. 

3.  Derivation  of  the  Nonlinear  Development  of  the  Instability 

From  observations  the  mean  square  density  fluctuation,  after  stabilization 

is  estimated  to  be  of  the  order  of  ' 3 x 10"^  (Booker,  1960). 

n 

Correspondingly  we  consider  a small  quantity  e,  say  10"^<e<10”^ , such  that 
we  may  write 


6 


5n  ' en 


(14) 


Also  the  experimental  observations  of  Kelley  and  Mozer  (1973)  show  that 

[E  1'|E  |.  This  means  that  the  primary  electric  field  is  of  the  same  order 
X y 

as  the  secondary  electric  field.  This  implies  that  the  auroral  electro jet 
has  a current  in  the  East-West  direction  which  is  about  the  same  order  as 
that  in  the  y-direction.  It  is  to  be  noted  that  the  current  in  the  y- 
direction  results  from  the  continuity  equation  involving  the  field-aligned 
currant.  Thus  we  may  conclude  that  the  electron  drift  velocity  ' l^y  ^ ’ 

Based  on  observations  in  the  auroral  ionosphere  we  now  introduce  the 
following  order  of  magnitude  approximations: 

V ' E n 

e e 


V.  - e V - 
i e e 


n ~ EV.  - E^Q 

i i e 


m - V.  - e^Q 
i e 


c - V 
s e 


(15) 


2itc 


1/2 

where  c^  » {kg(T^  + T^)/m^}  * ion  sound  speed, 

''l  ■ 

i i e 


06) 


where  T • T^  = T^  is  the  plasma  temperature,  and  — - 10”“*  in  order  of 
magnitude.  Using  the  order  of  magnitude  approximations  as  indicated  in 
(15)  we  find  from  equation  (13)  that 
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We  now  expand  the  governing  equation  (4)-(8),  and  ( i.  i i*-  »*•>«« 
the  relations  (15, -(18).  To  lowest  order,  it  is  found  :n*’ 
without  changing  form.  This  means  that  the  growth  r«t«»  j' 
and  also  the  dispersive  effects  are  small  compared  witt.  -am  « 
the  frequencies  and  proceed  on  larger  time  scales.  I*.  ?r 
order  in  e we  use  the  multiple  time  scales  expansion  *•- 
and  Krylov  (1947),  i.e.  let 


6n(t)  = n (t^;  t^^;  t2;-.-) 


where  tj^  » ^*’o’  ^2  ” 

Also  let  6n  » dn^°^  + e(Sn^^^+... 


Then, 


3 r /.  V - 3 r (o)  . . ,36n 


(o) 


jin 


From  now  on  the  time  variable  will  be  expressed  in  t« 


subscript  ^ will  be  dropped  although  it  is  alwavs  iapii* 


We  expand  equation  (7)  by  means  of  the  scaling  r«4.ati.m» 
compare  the  order  of  magnitude  of  each  Individual  term 


(e^,  £^,  e^,  E^,  e,  1) 


Hence,  we  find  that 


s . 6V  - 0 


-e 


■e  , 


S . 6v  - n -ly  i • (S  X 6V 
-e  e e z -e 


3 (o' 


From  eq.  (5)  with  j ■ e the  orders  of  magnitude  of  the  individual  terms  are 


(1.  1,  e’^  1) 


and  eq.  (5)  splits  into 


(s  • 5V  - 0 

3s  -e 


[_  + (V^  . s)  6n<°>  + n A (S  . ^v^d))  ^ ± [(s  • ]-  0 

\ 26) 


[(S  • 5V^^°^)6n^^^I  - 0 (27; 


On  the  other  hand  the  orders  of  magnitude  of  the  individual  terms  of  eq.  (8) 
with  J ■ 1(e)  are 


[l(c2),  c(c2),  e(e2),  1(e),  1(1),  1(1),  e(l)] 


where  the  orders  of  magnitude  of  the  electron  terms  are  given  in  parentheses. 
Hence,  eq.  (8)  becomes 


« • + v,6v/°h  - - + o i . S X 5V 
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and 


36V. 

r 


(1) 


36V 


(o) 


3t 


+ e- 


-i 


(o) 


(1) 


3t 


+ !!i  y g . <sv  _ n.  i_  . s X 6y/°> 


C ‘ _ (1) 


.(o) 


e 


-e 


3t 


+ Ji  i • § X 6V 
e z -e 


n 3t 

(1) 


1 z 


36n 


(1) 


Note  chat  Che  term  <- 


's  3s 


->  disappears  as 


(29) 


c 2 c 2 

<-2 -L  <6n^^^>  - 0 

A n 3s 


The  orders  of  magnitude  of  the  individual  terms  of  eq.  (5)  with  j ■ i are 


(1.  e,  1.  e) 


Hence,  eq.  (5)  with  j ■ i can  be  decomposed  into 


and 


36n^°^  . ' 3 (o),  „ 

—IT"  " 3l  • '^^i  > “ ° 


36n^^^  36n^°^ 

+ e- 


3t 


(30) 


3,  + 'ii  • h 


+ (§  • 6v/°^6n^°^)  - 0 

aS  “1 


(31) 


From  eqs.  (23),  (26),  (28),  and  (30)  we  obtain 


Eq.  (32)  represents  the  dispersion  relation  for  a density  wave  propagating 
in  a direction  perpendicular  to  the  ambient  magnetic  field.  We  now  let  Che 
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density  perturbation  be  of  the  form 


5 n^°^  ■ 6n^°^  exp  [1  (ujt  - ks)  ] 


where  w ■ + lu^  (ui^,  real) 


First  we  let  ■ o,  and  substitute  the  form  of  6n''  ' as  given  by  (33)  Into 
the  dispersion  relation  (32).  Ve  obtain  the  phase  velocity  of  the  wave 


if  - (L  • a)  (1  - 


Moreover  eq.  (32)  becomes: 


!£i)  + (V  . S)  -L  Snl-) 

n n/  3t  -e  ^ 3s  " 

e 1 


which  yields  the  solution 


5n^®^s;  t;  et;  ...)  - 5n^°^[s  - t(^^  • S)(l-  et;  ...] 

e 1 


which  shows  that  the  profile  of  the  density  wave  propagates  without  distortion 
on  the  "fast"  time  scale  under  the  condition  given  by  eq.  (34).  On  the  other 
hand,  the  condition  for  Instability  may  be  obtained  by  setting  < 0 In  the 
dispersion  relation  (32).  Thus  waves  with  phase  velocities  given  by 


s “f'  <i.  • - sir'’' 


will  be  unstable. 


Order  of  Magnitude  Estimate  of  the  Fluctuation. 
We  now  substitute  (36)  into  eq.  (30)  and  obtain 


i.e. 


1 3n 
n 


(o) 


3t 


- 7(V  • §)(1  - 

n 


Vi  3 . (o) 


3V  ^ (V  • S)  (1  - .^)"'6n 

-i  n “®  0. 

n g 


e"i 

(o) 


(38) 


This  relation  is  then  substituted  into  eq.  (28),  and  we  make  use  of  the 
threshold  relation  (3A).  Also  we  let 


s i + s i 

XX  y y 


5V  - 6V  i + 6V  £ 

-e  ex  X ey  y 


(39) 


(AO) 


Noting  that  § • dV  *0,  eq.  (28)  then  becomes 


(Al) 


Therefore 

<5n^°^6V  - (S  • £ ) 4 7^(V  • S)(l  - c-(o). 

ey  x'  n 0^^  -e 


or. 


<«„SV^y>  - <S  • • 9)(1  - <«»«"» 


The  orders  of  magnitude  of  the  terms  in  eq.  (2b)  with  J ~ i are 

(1,  e,  e^,  1,  e,  1). 

.th 


(A2) 


Hence,  to  zero  order  eq.  (2b)  becomes 


^Yi  \ 

^ — E - v.V^ 

3t  m^  - i-1 


(A3) 
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As  assumed  previously,  Che  average  quancicles  are  in  a state  of  equilibrium, 


so 


0 


Hence,  eq.  (43)  reduces  to 


V.  - — ^ E 
-1  - 


and, in  particular,  its  y-component  for  the  ions  becomes 


(44 


V » — - — E 

iy  y 

By  a similar  analysis  eq.  (2b)  with  j > 

E + - V X B - 
c -e 

In  component  form  we  thus  have 


e becomes  to  lowest  order 


0 


(45 


(46 


Hence,  V • s 


V 

ex 


V 

ey 


I "x 


(V  s + V s ) 

av  V ov  V 


(47 


E I 


X 


Combining  (48),  (42),  (43),  and  (13)  yields 


(§ 


i ) - E (s  • i )]  1 

y y X 

- n f — E , 

B y ’ 


V V. 

e i 


Q n,  , 

e ij 


-1 


<6n6n> 


(49) 


where  <6n6V^y>  has  been  neglected  compared  with  n 
Rearranging  the  above  equation,  we  obtain 


{E  (i  • i )(§  • i ) 
X X y 


n n/  nn  ^ 

e i 


(s 


V n. 


(1  - 


'^e'^iv~^  <6n6n>,~^ 


e i 


nn 


(50) 


i i iy 

The  secondary  electric  field  E « ^ decreases  due  to  the  fact  that  the 

j ye 

ion  flux  in  the  y-direction  decreases.  Hence,  according  to  the  first 

of  eq's  (47),  the  East-West  electron  flux  also  decreases  until  the  nonlinear 
density  perturbation  in  the  wave  reaches  a saturated  state. 


i.e. 
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(51) 


according  to  eq.  (34). 

Suppose  the  wave  propagates  in  the  East-West  direction.  Then, 
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and,  from  (50),  we  see  that. 


if  Ey  j*  0,  we  have 
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It  was  observed  by  Whalen  and  McDlarmid  (1972)  that  the  auroral  electron 
precipitation  Is  practically  field-aligned  as  the  rocket  payload  passed  through 
the  Northern  edge  of  a visible  auroral  display.  Based  on  this  observation  we 
can  determine  the  magnitude  of  the  electron  precipitation  by  calculating  the 
field-aligned  current.  The  latter  is  related  to  the  horizontal  current  by 
the  current  continuity  equation  7 • j > 0.  Thus  we  have 


where 


a E + a„  E 

p y H X 


Op  and  Ojj  are  the  height  Integrated  Pedersen  and  Hall  conductivity,  respectively, 
and  j , represents  the  field-aligned  current  density. 


Substituting  eq.  (50)  into  (53),  we  obtain 
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If  now,  for  convenience,  we  assume  that  this  current  Is  uniform  over  the 
arc  width,  then  we  may  determine  the  field-aligned  current  density  . 
Consequently  the  electron  precipitation  may  be  calculated. 


5.  Conclusions 


Based  on  the  Farley-Buneman  two-stream  instability  (Buneman,  1963; 

Farley,  1963)  we  have  developed  a nonlinear  analysis  of  the  auroral  electrojet. 
The  electron  temperature  is  assumed  to  be  equal  to  the  ion  temperature  in 
the  electrojet.  The  results  of  our  analysis  indicate  the  following  points 
of  Interest. 

1)  The  nonlinear  analysis  of  our  two-fluid  model  is  able  to  predict  the 
order  of  magnitude  of  the  irregularities  of  the  auroral  electrojet. 
Obviousl>  this  is  outside  the  scope  of  a linear  theory.  If  we  use  the 
results  of  Kelley  and  Mozer's  (1973)  observation  that  the  primary 
electric  field  is  approximately  equal  to  the  secondary  electric 
field  Ey  (in  the  South-North  direction) , then 
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This  fluctuation  level  is  of  the  same  order  of  magnitude  as  that  observed 
by  Booker  (1960) . 

2)  The  nonlinear  theory  shows  that  the  electron  drift  velocity  for  instability 
is  consistently  larger  than  the  phase  velocity  of  the  unstable  wave.  The 
latter  is  slightly  larger  than,  or  equal  to,  the  ion-acoustic  velocity. 

This  result  is  in  contrast  to  the  linear  theory  of  the  Farley-Buneman 
instability,  where  the  drift  velocity  is  shown  to  be  larger  chan  the  ion 
thermal  velocity. 

3)  From  eq.  (50)  we  note  Chat  when  the  direction  of  the  wave  propagation 
la  exactly  in  the  y-direction,  then  E^  is  reduced  to  zero.  This  not 
only  does  not  agree  with  Kelley  and  Mozer’s  observations  which  indicates 
a non-zero  E^  field,  but  also  contradict  our  assumption  of  the  Blrkeland 
current  model.  Therefore,  we  conclude  that  Che  wave  propagation  vector 
should  have  a small,  but  non-zero,  component  in  Che  x-dlrectlon. 
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